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Detection of superimposed periodic signals using wavelets 
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1 INTRODUCTION 



ABSTRACT 

In this paper we present a wavelet based algorithm that is able to detect superimposed 
periodic signals in data with low signal-noise ratio. In this context, the results given by 
classical period determination methods highly depend on the intrinsic characteristics of 
each periodic signal, like amplitude or profile. It is then difficult to detect the different 
periods present in the data set. The results given by the wavelet based method for 
period determination we present here are independent of the characteristics of the 
signals. 

Key words: methods: data analysis - stars: variables: other - X-rays: general. 
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of astro nomy, and, often, simple Fourier based methods or 

epoch folding methods are able to detect the periods that 
are present in the data sets. However, this turns out to be 
difficult when the periodic signals are not sinusoidal or when 
the signal-noise ratio is very low. 

Specifically, in high energy astrophysics two problems 
are responsible for reduced signal-noise ratios in the data: 
the relatively small number of high energy photons when 
compared to the number of photons available at other wave- 
lengths, and the low efficiency in photon counting of high en- 
ergy detectors. Hence, the astronomical data supplied by X- 
ray telescopes, especially flux monitorings, often suffer from 
poor statistics. This poses many problems when processing 
and analysing the data in order to determine the flux, detect 
periodic signals, etc. In addition, some X-ray sources present 
several periods, which arise from different physical processes 
such as pulses, eclipses in two-body systems or occultation 
by a precessing accretion disc. Some examples of this kind 
of sources are SMC X-l ( |Wojdowski et al 199&1), LMC X-4 
and LS I +65°010 flCorbet, Finley & 



1989; 



Dreher 1987| and press fc Rybicky 1989| ) 



Roberts, Lehar & 



1991) 



( Levine et al 
Peele 1999| ).~ 

Classical period determination methods are based ei- 
ther on epoch folding techniques or on Fourier decompo- 
sition analysis. The first group of methods are based on 
the analysis of the dispersion of the diferent light curves 
produced by foldi ng the data over a range of trial periods 
(see, for example, Lafler fc Kinman 1965| ; Jurkevich 1971 



and Btellingwcrf 1978| ). The second type of methods use the 



Fourier transform in combination with deconvolution tech- 
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Epoch folding methods and Fourier based methods work 
well when applied to data sets that present a unique period. 
However, when the analysed data set contains several peri- 
odic signals, the behavior of these algorithms highly depends 
on intrinsic signal characteristics. Every individual periodic 
signal may present different spectral behavior (amplitude 
and profile) and depending on them, the algorithm will de- 
tect some kinds of signal better than others. Fourier based 
techniques will, in general, successfully detect two combined 
sinusoidal signals, even if one of them has a low signal-noise 
ratio. However, they have difficulties in detecting the non- 
sinusoidal signals that might be present in a data set (as in 
the astronomical case of pulsed emission superimposed on 
an orbital period). On the contrary, epoch folding methods 
are well suited for detecting non-sinusoidal signals but they 
tend to fail when two or more periods are present in the 
data set, especially when the signal-noise ratio is low. 

In general, we can say that the greater the difference 
between the spectral characteristics of each signal, the more 
difficult it is to detect each period with classical methods, 
and the less statistically significant are the results given by 
them. Therefore, a new algorithm becomes immediately nec- 
essary for the detection of superimposed periodic signals. 

In this paper, we show how the methodology of wavelet 
theory is very well suited to this problem, since it is com- 
pletely oriented to decompose functions into their several 
spectral characteristics. This allows us to isolate every signal 
present in our data and to analyse them separately, avoid- 
ing their mutual influences. In Section 2 we outline some 
concepts in wavelet theory relevant to the stated problem. 
In Section 3 we propose an algorithm to detect each of the 
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periodic signals present in a data set by combining wavelet 
decomposition with classical period determination methods. 
In Sections 4 and 5 we present some examples of synthetic 
data we used to test the algorithm and the results we ob- 
tained. We summarise our conclusions in Section 6. 



2 OUTLINE OF THE WAVELET TRANSFORM 

Multiresolution analysis based on the wavelet theory intro- 
duces the concept of details between successive levels of scale 



or resolutions (Chui 1992 


; Daubechies 1992 


; |Meyer 1993 


Young 1993; 


Kaiser 1994; 


Vetterli & Kovacevic 1995|). 



Wavelet decomposition is i ncreasingly being used for 
astronomical signal processing ( [gtarck fc Murtagh 1994 



Starck, Murtagh fc Bijaoui 1998) and rem otely sensed data 



(yocky 1991; Datcu, Luca fc Seidel 199Q. The method is 



based on the decomposition of the signal into multiple chan- 
nels based on their local frequency content. The wavelet 
transform is an intermediate representation between the 
Fourier and the temporal one. In the Fourier transform we 
know the global frequency content of our signal, but we have 
no information about the temporal localization of these fre- 
quencies. However, the wavelet transform gives us an idea 
of both the local frequency content and the temporal distri- 
bution of these frequencies. Since in Fourier space the basis 
functions are sinusoidal, they extend along all space and 
do not have temporal concentration. However, wavelets are 
concentrated around a central point, thus they have a high 
degree of temporal localization. 



2.1 Some theory 

Next, we will briefly present an outline of the wavelet trans- 
form, relying on the multiresolution signal representation 
concept. 

Given a signal f(t), we construct a sequence F m (f(t)) 
of approximations of f(t). Each F m (f(t)) is specific for the 
representation of the signal at a given scale (resolution). 
Fm(f(t)) represents the projection of the signal f(t) from 
the signal space S onto subspace S m ■ In this representation, 
Fm(f(t)) is the 'closest' approximation of f(t) with resolu- 
tion 2 m . 

The differences between two consecutive scales m and 
m + 1 are the multiresolution wavelet planes or 'detail' signal 
at resolution 2 m : 



»m(/(i)) = F m (f(t)) - F m+1 (f(t)) . 
This detail signal can be expressed as: 

Wm(/(t)) =£V m ,l(/)V>m,l(t) , 



(1) 



(2) 



where coefficients W m ,i(f) are given by the direct wavelet 
transform of the signal f(t): 



W m ,l(f) 



f(t)ll> m ,t(t) dt 



(3) 



The coefficients W m ,i(f) are called wavelet coefficients 
of f(t). Such coefficients correspond to the fluctuations of 
the signal f(t) near the point I at resolution level m. Thus, 
the wavelet transform (H) represents the expansion of signal 



f(t) in the set of basis functions ip m ,l(i). These basis func- 
tions are scaled and translated versions of a general function 
ip(t) called Mother Wavelet. To construct the basis func- 
tions ip m ,i{t), the Mother Wavelet is dilated and translated 
according to parameters m and I as follows: 

VwW = 2 m/ V(2 m i-0 . (4) 

Therefore, all the basis functions ip m ,i(t) have the same pro- 
file, that is, the Mother Wavele t profile. Using (Iff ) we obtain 
an orthonormal wavelet basis ( Daubechies 1992] ). 

The inverse discrete wavelet transform is given by the 
reconstruction formula: 



/(*) = £S w wtfw™.*(*) 



(•>) 



In summary, the wavelet transform describes at each 
resolution step the difference between the previous and the 
current resolution representation. By iterating the process 
from the highest to the lowest available resolution level we 
obtain a pyramidal representation of the signal. This usually 
includes a decimation process, i.e., in each iteration 1 out of 
2 points is discarded, which implies that the number of data 
points at lower frequencies is highly reduced. 

In the Fourier transform, the noise contribution is 
spread along all frequencies. In contrast, one of the inter- 
esting properties of the wavelet transform (which is also fre- 
quency based) is that noise contribution is more important 
on the higher frequency wavelet planes. 



2.2 The 'a trous' algorithm 

In order to ob tain a discrete wavelet de composition for sig- 
nals, we follow Starck & Murtagh (1994) and we use the dis- 
crete wavelet transform algorithm known as 'a trous' ('with 
holes') to decompose the signal into wavelet planes. Given 
a signal p we construct the sequence of approximations: 



Fi(p)=Pi, F 2 (pi)=p2, F 3 (p 2 )=p 3 ,- 



(6) 



To construct the sequence, this algorithm performs suc- 
cessive convolutions with a filter obtained from an auxiliary 
function named scaling function. We use a scaling function 
which has a B3 cubic spline profile. The use of a B3 cubic 
spline leads to a convolution with a mask of 5 elements, all 
elements being scaled up by 16: (1,4,6,4,1). 

As stated above, the wavelet planes are computed as 
the differences between two consecutive approximations pi-i 
and pi. Letting Wi = p%-\ — pi (i = 1, ■ ■ ■ , n), in which 
Po = P, we can write the reconstruction formula: 



^ Wi + p r 



(7) 



In this representation, the signals pi (i — 0, • ■ • , n) are 
versions of the original signal p at increasing scales (decreas- 
ing resolution levels), Wi (i = 1, • ■ • , n) are the multiresolu- 
tion wavelet planes and p r is a residual signal (in fact n = r, 
but we explicitly substitute n by r to clearly express the 
concept of residual). In our case, we are using a dyadic de- 
composition scheme. Thus, the original signal po has double 
resolution than pi, the signal pi double resolution compared 
to P2, and so on. If the resolution of signal po is, for example, 
10 At (being At the sampling of the data), the resolution of 
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pi would be 20 At, the resolution of P2 would be 40 At, etc. 
All these pi (i = 0, • • • , n) signals have the same number 
of data points, in contrast to some more usual wavelet de- 



detections present in PDM and CLEAN and to detect the true 
periods present in our original data. 



We can use ot her wavelet based algorithms (Szatmary, 



composition algorithms which reduce the number of points 



Vinko & Gal 1994), which are based on approximations of 



by a factor 2 when going through increasing scales (process 
known as decimation step). Later we will see why we are 
interested in maintaining the number of points. 



3 PROPOSED ALGORITHM 

We propose to apply the wavelet decomposition using the 
'a trous' algorithm to solve our initial stated problem: to 
isolate each of the periodic signals contained in a set of data 
and study them separately. 

Hence, the algorithm we propose is as follows: 

(i) Choose a value for n and decompose the original signal 
p into its wavelet planes Wi (i = 0, • ■ • ,n). 

(ii) Detect periods in each of the n obtained wavelet plane 

VJi. 

Any method can be used to detect the periods present 
in each Wi. In ou r tests, we used Phase Dispersi on Min- 
imization ( pdm) ( Stellingwerf 1978) and CLEAN (Roberts 



et al. 1987) methods. PDM belongs to the group of epoch 



folding methods mentioned in the introduction. Therefore it 
folds the data points over a set of trial periods and analyses 
the dispersion of every light curve in the phase space. To 
do that, the data points are grouped in phase bins and the 
most probable period is then found by minimizing the sum of 
the bin variances. On the other hand, CLEAN belongs to the 
group of Fourier transform decomposition methods. It works 
in the frequency domain and it is based on the fact that the 
spectrum obtained from the Fourier transformation of the 
data set contains artefacts caused by the imcompleteness 
of the sampling (time gaps) and the finite time span of it. 
CLEAN uses a non-linear deconvolution algorithm to clean 
(hence its name) these artefacts from the original (dirty) 
spectrum. This is done by interatively substracting from the 
dirty spectrum the expected response of a signal composed 
by a unique harmonical function (clean component). From 
the set of clean components, a clean spectrum is obtained 
which is at the end corrected to preserve the noise level and 
the frequency resolution. 

PDM is well suited to identify signals with any profile, 
not just a sinusoidal one. For instance, if we have a burst-like 
periodic signal, it is more successfully detected by the PDM 
algorithm than by CLEAN. On the other hand, if we knew 
that we are trying to estimate the period of a sinusoidal 
signal, it would be better to use CLEAN. 

The algorithm we present here could be seen as an im- 
provement over the usual period detection algorithms, be- 
cause prior to using them we decompose our signal into its 
Wi wavelet planes, and afterwards we apply these algorithms 
to each Wi. 

In PDM and CLEAN methods there is only one data set 
to study and to estimate its period. But in the wavelet based 
method there are several data sets to study (several wavelet 
planes), and this improves the detection probabilities. We 
can search for a certain period in several wavelet planes, 
which helps to discard some of the usual marginal artefact 



a continuous wavelet transform and the study of wavelet 
space coefficients. However, these algorithms present a non 
direct inverse wavelet transform: the search for periodicities 
is based on the fit between the wavelet base function profile 
and the signal one, and therefore on the values of the wavelet 
transform coefficients, which highly depend on the wavelet 
base used. For example, if we took a signal constructed by a 
Gaussian with 1 day duration which repeats every 30 days, 
it would be very difficult to find a wavelet base that fits 
well to this sparse signal. On the contrary, if we take a suit- 
able wavelet base and a decomposition scheme which allows 
us to perform easily its inverse transform, this base will fit 
the profile of every Gaussian. Hence, performing the inverse 
transform of every wavelet scale, we can work in the tempo- 
ral representation to find directly the periodicities. This is 
the case of the wavelet decomposition algorithm we use in 
the method we present here. 

Another reason to use the 'a trous' algorithm is that the 
sampling of the wavelet planes is the same as the original 
data. This allows us to work directly in the temporal space 
with the frequency content of the corresponding wavelet 
plane. In other wavelet decomposition algorithms, we are 
usually forced to work on the decimated wavelet space. 

Hereafter, and for notational convenience, the wavelet 
based PDM and CLEAN methods will be called WPDM and 
WCLEAN respectively. 



4 SIMULATED DATA 

In order to check the benefit of applying WPDM versus PDM, 
or WCLEAN versus CLEAN, we generated several sets of sim- 
ulated data containing two superimposed periodic signals. 
Each data set is composed of a high amplitude primary si- 
nusoidal function and a secondary low amplitude one. For 
the latter we used two kinds of functions: sinusoidal and 
Gaussian. The first are intended to simulate variable sources 
with pure sinusoidal intensity profile (like precession of ac- 
cretion discs), and the second burst-like events (like pulses 
or eclipses). Finally, we added a white-Gaussian noise to this 
combination of signals. We increased the value of the noise 
standard deviation, a, up to the value where detection of pe- 
riodic signals became statistically insignificant in both the 
classical and the wavelet based methods. 

Here we list the characteristics of the signals: 

(i) Each signal is generated as an evenly spaced data set 
composed of 1000 points, separated by a unit between them. 

(ii) The high amplitude sinusoidal function has an am- 
plitude equal to 1 and a period of 108.5 units. The period 
selected should satisfy two conditions: not be an exact divi- 
sor of the number of points in the data set, and allow the 
presence of more than 5 complete periods in the simulated 
signal to guarantee the possibility of detection. 

(iii) The amplitudes of the low amplitude periodic func- 
tions (sinusoidal or Gaussian) are 0.1 and 0.5. 

(iv) The periods used for the secondary functions are 
13.13 and 23.11 units. Others could be used, satisfying the 



4 X. Otazu et al. 



SIMULATED DATA 



3.00 




400 600 
Time units (At) 



1000 
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CLEAN 
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Trial Period (At) 



50 100 
Trial Period (At) 



150 



Figure 1. Top: simulated data of a sinusoid with 23.11 units 
period, amplitude=0.1 and cr noise = 0, superimposed on the si- 
nusoid with 108.5 units period and amplitude=1.0. Bottom: PDM 
and CLEAN outputs of this data set. 



condition that they are not harmonics of the high amplitude 
sinusoidal function period. 

(v) In the case of the Gaussian signal, we have used two 
values for the Full Width Half Maximum (FWHM), corre- 
sponding to 2 and 6 units, respectively. 

In the first three columns of Table [l] we present the 
parameters used to generate each simulated data set in 
the case of sine+sine periodic signals. In Table ^| the first 
four columns show the parameters used to generate the 
sine+Gaussian signals. 



5 RESULTS 

The simultaneous use of two independent methods, such as 
CLEAN and PDM, is usually applied to discriminate false pe- 
riod detections from the true ones. A similar procedure can 
be used with each one of the wavelet planes in the WPDM 
and WCLEAN methods. Therefore, when comparing the be- 
havior of these different methods, we have to compare the 
usual PDM-CLEAN method combination for period estima- 
tion prior to the new WPDM- WCLEAN combination. 

The primary period (108.5 units) is always detected b" 
all methods, and does not appear in Table 1 and Table 
In these tables, the four last columns show the detected 
low amplitude periods for each data set using PDM, CLEAN, 
WPDM and WCLEAN methods, respectively. A dash is shown 
when a period is not detected, and a question mark when the 
detection is difficult or doubtful. When a period is indicated 
in the wavelet based methods, we also show in parentheses 
the wavelet planes where it is detected. 



5.1 Sine+sine case 

To illustrate the performances of WPDM and WCLEAN meth- 
ods, we show in Fig. |l| a simulated data set which contains 
two pure sinusoidal functions without noise, and both their 
PDM and CLEAN outputs (in PDM we look for minima with 
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Figure 2. WPDM output for each wavelet plane of the data at 
the top of Fig. [l| 




Trial Period (At) 

Figure 3. WCLEAN output for each wavelet plane of the data at 
the top of Fig. [l| 



subharmonics, while in CLEAN we look for maxima). The 
108.5 units period is detected by both methods. In the case 
of PDM we have zoomed into the range of short trial periods 
and dispersion to try to detect the 23.11 units minimum. It 
is clear that even in the case of signals without noise, PDM 
method is unable to find the secondary signal (period=23.11 
units). In contrast, CLEAN is able to find it, although with a 
low normalized power (less than 0.1, which is the ratio be- 
tween the amplitudes of both periodic signals). The output 
of WPDM is shown in Fig. [] and the one of WCLEAN in Fig. [| 
As one would expect, the secondary period is only seen in 
the lower wavelet planes (wi to W4), while the primary long- 
term period of 108.5 units is better detected when working 
with higher wavelet planes (w± to w§). We have worked un- 
til we, because in this plane the 108.5 units period is clearly 
detected, and there is no reason to continue with higher val- 
ues of n. This will be the case in all the analysed data we 
present here. 

Some remarks can be made after inspection of Table |l| 

(i) There is a better performance of CLEAN relative to 
PDM, because the latter has problems when dealing with 
superimposed low amplitude periodic signals. Only when the 
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Table 1. Periods detected in the sine+sine data sets. The first three columns are the 
parameters used for the simulated data. A dash is shown when a period is not detected, 
and a question mark when the detection is difficult or doubtful. When a period is indicated 
in the wavelet based methods, we also show in parentheses the wavelet planes where it is 
detected. 



Period Amplitude cr n 



PDM 



CLEAN WPDM 



WCLEAN 



13.13 



0.1 



0.5 



23.11 



0.1 



0.5 



0.0 
0.1 
0.2 
0.4 
0.6 
0.8 
1.0 
1.2 
1.4 
1.6 
0.0 
0.5 
1.0 
2.0 
3.0 
4.0 
5.0 
6.0 
8.0 
0.0 
0.1 
0.2 
0.4 
0.6 
0.8 
1.0 
1.2 
0.0 
0.5 
1.0 
2.0 
3.0 
4.0 
5.0 
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23.11 
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23.12? 



23.11 
23.11 
23.12 
23.11 
23.11? 



13.13 
13.13 
13.11 
13.09 
13.09 
13.08 
13.08 
13.07 
13.07 

13.13 
13.13 
13.13 
13.13 
13.13 
13.13 
13.13 
13.13 

23.11 
23.11 
23.13 
23.11 
23.17 
23.13 
23.13 

23.11 
23.11 
23.12 
23.11 
23.17 
23.13 



1,2,3) 

2,3) 

2,3) 

2,3) 

3) 

3) 

3) 

3?) 

3?) 

1,2,3) 

2,3) 

2,3) 

2,3) 

2?,3) 

3) 

3) 

3?) 

1,2,3,4) 

2,3,4) 

3,4) 

3,4) 

3,4) 

3?,4) 

4?) 

1,2,3,4) 

2,3,4) 

2?, 3,4) 

3,4) 

3?,4) 

4?) 



13.13 
13.13 
13.11 
13.09 
13.09 
13.08 
13.08 
13.07 
13.07 

13.13 
13.13 
13.13 
13.13 
13.13 
13.13 
13.13 
13.13 

23.11 
23.11 
23.13 
23.11 
23.17 
23.13 
23.11 

23.11 
23.11 
23.13 
23.11 
23.17 
23.13 



1,2,3) 

2,3) 

2,3) 

2,3) 

2,3) 

2,3) 

2,3) 

2?,3?) 

2?,3?) 

1,2,3) 

2,3) 

2,3) 

2,3) 

2,3) 

2,3) 

2?,3) 

3?) 

1,2,3,4) 

2,3,4) 

3,4) 

3,4) 

3,4) 

3?,4) 

4?) 

1,2,3,4) 

1,2,3,4) 

2?,3,4) 

3,4) 

3?,4) 

4?) 



amplitude of the secondary signal is close enough to the 
primary one (amplitude=0.5) is PDM able to detect it. 

(ii) In all methods, with high noise-signal ratios the de- 
tected periods are slightly different from the simulated ones. 

(iii) WPDM and WCLEAN perform always better than 
PDM and CLEAN methods. When CLEAN marginally detects 
the secondary period, WPDM and WCLEAN have no prob- 
lems to detect it, and they work properly even with higher 
noise. In the WPDM case, the results are always much better 
than with PDM. 

(iv) As the noise increases, the detection starts to fail 
in the lower wavelet planes (higher frequencies), and only 
the higher ones (lower frequencies) are noise-free enough to 
allow period detection. 

(v) The maximum ratio between the noise and the signal 
amplitude to detect a period with WPDM and WCLEAN, is 
nearly constant for a given period, being around 12 for the 
13.13 units period and around 8 for the 23.11 units period. 
The reason for this difference is that there are 76 complete 
periods in the case of 13.13 units period, and only 43 com- 
plete periods in the case of 23.11 units period. 
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Figure 4. Wavelet decomposition of the simulated data which 
contains a Gaussian with 13.13 units period, a FWHM of 2.0 
units, amplitude=0.1 and cr no i ae = 0, superimposed on the sinu- 
soid with 108.5 units period and amplitudc=1.0. 
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Table 2. Periods detected in the sine+Gaussian data sets. The first four columns are the parameters 
used for the simulated data. A dash is shown when a period is not detected, and a question mark when 
the detection is difficult or doubtful. When a period is indicated in the wavelet based methods, we also 
show in parentheses the wavelet planes where it is detected. 
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5.2 Sine+Gaussian case 

For illustrative purposes, we show in Fig. ^ the wavelet de- 
composition of a sine+Gaussian signal. In the lower wavelet 
planes (planes u>i and W2) we have isolated the Gaussian 
contribution. In the w$ we have a mixed contribution from 
both periods, but in the higher planes we only have the 
longest period contribution. 

We must note that the use of two different FWHM for 
the Gaussian, combined with two different periods (13.13 
and 23.11 units), gives 4 different profiles. Hence, the phase 
duration of the burst-like event ranges from very low to rel- 
atively high values in the following order: FWHM=2 and 
period=23.11, FWHM=2 and period=13.13, FWHM=6 and 
period=23.11, and finally FWHM=6 and period=13.13. 

In view of the results displayed in Table ^ we can make 
the following comments: 

(i) Again, there is a better performance of CLEAN over 
PDM, for the reason explained above. However, we must note 
that when the FWHM is only 2 units, PDM never detects 
the secondary period. Only with FWHM=6 units and an 
amplitude of 0.5, can PDM detect the low amplitude periodic 
signals. 



(ii) In all methods, with high noise-signal ratios the de- 
tected periods are slightly different from the simulated ones. 

(iii) WPDM and WCLEAN perform always better than 
PDM and CLEAN methods. When CLEAN marginally detects 
the secondary period, WPDM and WCLEAN have no prob- 
lems to detect it, and they work properly even with higher 
noise. In the WPDM case, the results are always much better 
than with PDM. 

(iv) As the noise increases, the detection starts to fail 
in the lower wavelet planes (higher frequencies), and only 
the higher ones (lower frequencies) are noise-free enough to 
allow period detection. 

(v) In all cases with amplitude=0.1, the WPDM perfor- 
mance is very similar to WCLEAN. In the amplitude=0.5 
cases, wpdm is always better than WCLEAN because the 
signal is not sinusoidal and the amplitude is high enough to 
allow a good detection. 

(vi) For a given amplitude of the Gaussian signal, the 
maximum noise-signal ratio achieved with WPDM and 
WCLEAN increases with the phase duration of the FWHM. 

Finally, and for illustrative purposes, we show in Fig. ^ 
and Fig. H two simulated data sets generated with the fol- 
lowing common parameters: 13.13 units period, FWHM=6.0 
units and amplitude=0.5. The only difference is that in 
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Figure 5. Top: simulated data of a Gaussian with 13.13 units 
period, a FWHM of 6.0 units, amplitude=0.5 and cr no i 3e = 0.5, 
superimposed on the sinusoid with 108.5 units period and ampli- 
tude^. 0. Bottom: PDM and CLEAN outputs of this data set. 
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Figure 8. Top: simulated data of a Gaussian with 13.13 units 
period, a FWHM of 6.0 units, amplitude=0.5 and cr no i ae = 1.5, 
superimposed on the sinusoid with 108.5 units period and ampli- 
tude^. 0. Bottom: PDM and CLEAN outputs of this data set. 
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Figure 6. WPDM output for each wavelet plane of the data at 
the top of Fig. || 



Figure 9. WPDM output for each wavelet plane of the data at 
the top of Fig. || 
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Figure 7. WCLEAN output for each wavelet plane of the data at 
the top of Fig. || 



Figure 10. WCLEAN output for each wavelet plane of the data 
at the top of Fig. ^. 
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Fig. |[ a noise = 0.5, while in Fig. |, a noise = 1.5. The out- 
puts of PDM, CLEAN, WPDM and WCLEAN are shown from 
Fig. | to Fig. [o| 

As can be seen in the figures, in the case cr n0 i se = 0.5, 
the 13.13 units period is clearly detected by all methods. In 
WPDM and WCLEAN the period is detected in the wavelet 
planes u>2 and 103. These detections are much better than 
those of PDM and CLEAN. 

On the other hand, in the a n0 i se = 1.5 case, the period 
is not detected by PDM nor by CLEAN. WCLEAN also fails 
to detect it: the maximum in the wavelet plane 103 is not 
indicative of a period detection and, moreover, its position 
is in the 15.5 units trial period. The only clear detection 
with a subharmonic is in W3 of WPDM. 



6 CONCLUSIONS 

When dealing with two superimposed periodic signals, the 
performance of classical methods, like PDM and CLEAN, de- 
pends on the characteristics of these signals. These methods 
often fail to detect simultaneously both periods, even with 
high or moderate signal-noise ratios. 

With the wavelet based methods we present here, 
WPDM and WCLEAN, it is possible to reach lower values 
of signal-noise ratios and still detect both periodic signals. 
This means that the wavelet based methods are less affected 
by noise than PDM and CLEAN, and they allow us to detect 
periodic signals with noisier data. 

Another advantage of the proposed methods is the si- 
multaneous detection of periods in several wavelet planes. 
If a period is marginally detected in more than one wavelet 
plane, we are probably detecting a true periodic signal. This 
allows us to improve the detection confidence. 

These facts prove the major ability of WPDM and 
WCLEAN over PDM and CLEAN to deal with noisy signals 
and to detect superimposed periods even with low signal- 
noise ratios. 
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